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The cosmological analysis of X-ray cluster surveys: 

II- Application of the CR-HR method to the XMM archive 
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ABSTRACT 

We have processed 2774 high-galactic observations from the XMM archive (as of May 
2010) and extracted a serendipitous catalogue of some 850 clusters of galaxies based on 
purely X-ray criteria, following the methodology developed for the XMM-LSS survey. 
Restricting the sample to the highest signal-to- noise objects (347 clusters), we perform 
a cosmological analysis using the X-ray information only. The analysis consists in 
the modelling of the observed colour-magnitude (CR-HR) diagram constructed from 
cluster instrumental count-rates measured in the [0.5-2], [1-2] and [0.5-1] keV bands. 
A MCMC procedure simultaneously fits the cosmological parameters, the evolution of 
the cluster scaling laws and the selection effects. 

Our results are consistent with the eg and fl m values obtained by WMAP-5 and 
point toward a negative evolution of the cluster scaling relations with respect to the 
self-similar expectation. We are further able to constrain the cluster fractional radius 
Xc.o — r c/ R500C to x c _o = 0.24 ± 0.04. This study stresses again the critical role 
of selection effects in deriving cluster scaling relations, even in the local universe. 
Finally, we show that CR-HR method applied to the eRosita all-sky survey - provided 
that cluster photometric redshifts are available - will enable the determination of the 
equation of state of the dark energy at the level of the DETF stage IV predictions; 
simultaneously, the evolution of the cluster scaling-relations will be unambiguously 
determined. 

The XMM CLuster Archive Super Survey (XCLASS) serendipitous cluster catalogue 
is available online at: http://xmm-lss.m2p3. fr:8080/14sdb/. 

Key words: Cosmology: observations - Galaxies: clusters: general - X-rays: galaxies: 
clusters - Methods: observational - Catalogues 



1 INTRODUCTION 

Clusters of galaxies, the most massive bound objects in the 
universe, are the direct products of the growth of cosmic 
structures. Using cluster samples in cosmological analyses 
requires not only to span a large range of redshifts and 
masses. It is also mandatory to precisely understand how 
those objects have been selected and how the selection is re- 
lated to the cosmological distribution of dark matter haloes, 
the only quantity handled by the theory. Indeed, under- 
standing selection processes turned out to be one of the 
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main challenges of today's cluster cosmology and is inti- 
mately related to our ability to adequately determine clus- 
ter scaling-relations along with the associated dispersion. In 
this respect, X-ray surveys are potentially extremely pow- 
erful, given that the X-ray properties of the cluster pop- 
ulation can be derived from ab initio models. Substantial 
efforts have been devoted to assemble statistically signifi- 
cant cluster sam ples with the past generation of X-ray ob- 



servatories: (e.g. IScharf et al.lll997l: Vikhlinin et al.lll998a: 


Jones & Forman 1999 


: Bohrineer et all 200d: Borsani et alj 


200 ll: Bohringer et al. 


120041: Burenin et al.l2007). Ten vears 



ago, XMM opened a new era in cluster surveys, allow- 
ing us to access and to characterise clusters, much fainter 
than enabled by e.g., ROSAT. Nowadays, cluster serendip- 
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itous searches in the XMM archive arouse a growing inter- 
est thanks to the wealth of pointed observation s, publicly 
available : The XCS survey (jRomer et alJl200ll ). launched 
a decade ago is n ow delivering its first X-ray selected cata- 
logue s of clusters (|Llovd-Davies et al.ll20lj ; Mehrtens et al.l 
l201ll ). Se veral other projects are b eing conducted such as 
SExclass l|Kolokotronis et alJ 2006]) and combined searches 
with Chandr a archival data ([Peterson et al, 20091 ) or SDSS 
optical data (jTakev. Schwope. fc Lamerll201ll )! 

This paper is the second of a series describing a 
novel approach to the cosmological interpretation of clus- 
ter number counts in X-ray cluster surveys. The first pa- 
per l|Clerc et al.ll201ll . hereafter paper I) investigated the 
constraining power of a method based on the analysis of 
instrumental X-ray observables, namely a count-rate (CR) 
and a hardness-ratio (HR) in well defined X-ray bands. The 
combination of the two quantities was shown to reliably 
describe the surveyed cluster population: the correspond- 
ing CR-HR statistical distribution, which is analogous to 
a colour-magnitude diagram, can be fully predicted by an 
ab initio modeling involving the cosmology, cluster scaling- 
relations, the survey selection effects along with the XMM 
instrumental response. In this paper, we present an inde- 
pendent analysis of 2774 high-galactic latitude observations 
from the XMM archive having effective exposure times of 
10 and 20 ks. Followin g a selection proce dure adapted from 
the XMM-LSS s urvey flPierre et alj|2007t >. we detect 845 CI 



galaxy clusters l|Pacaud et al.l 20071 ) , hence constituting the 
X-CLASS catalogue (XMM Cluster Archive Super Survey). 
We apply the CR-HR method to a subsample of 347 clusters 
selected for their high signal-to-noise ratio over an effective 
geometrical area of 90 sq.deg. We devote special care to the 
count-rate measurements and to the derivation of the selec- 
tion function associated to this heterogeneous archival data 
The structure of this paper is as follows. We first present 
the steps leading to the creation of the X-CLASS cluster 
sample (Sect.0. In a second part (Sect. [3} we describe the 
derivation of the survey selection function and how we ac- 
count for the presence of pointed clusters in the sample. In 
Sect. U and [5] we present the results of the analysis of the 
CR-HR distribution. We discuss our results in Sect. [6] and 
present some cosmological predictions for the eRosita all-sky 
survey. Summary and conclusion are gathered a in Sect. [7] 



2 THE X-CLASS CATALOGUE 

The XMM CLuster Archive Super Survey (X-CLASS) is 
based on the analysis of archival observations from the 
XMM-Newton observatory. It is intended to provide a sam- 
ple of several hundreds of clusters suitable for cosmological 
studies. In this section, we describe the selection of the orig- 
inal data and the adopted methodology for detecting and 
characterizing galaxy clusters. 



2.1 Pre-selection of the XMM archival data 

As of May, 26th 2010, 7716 individual obse rvations were 
listed in the XMM Science Archive system (|Arviset et al.l 
2003). Out of these, we retrieved 2774 observations selected 
as follows: 



All Observations 
XCLASS processed 
XCLASS processed & filtered 




10 100 
Observation exposure time Texp (ks) 



Figure 1. Cumulative exposure distribution of the XMM archival 
observations. Dashed line: nominal exposure time for all 7716 
observations. Dot-dashed line: nominal exposure time for the 2774 
retrieved observations . Plain line: effective "clean" exposure time 
after processing and background flare removal. 
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Figure 2. Distribution of hydrogen column densities, (dashed 
line:) for the 7716 XMM archival observations, (plain line:) for 
the 2774 pre-processed observations. 



(i) pointing center at high galactic latitude (|6| ^ 20deg) 
to minimize the effect of galactic absorption, 

(ii) total exposure time greater than 5 ks, 

(iii) all three detectors (MOS1, MOS2 and PN) in imag- 
ing mode and at least one of them in Full Frame mode, 

(iv) pointing center not closer than 5 deg to the Small 
and Large Magellanic Clouds, and 2 deg away from M31, 

(v) public data (as of May, 26th 2010) 

This selection contains in particu lar 92 observation s 
from the 10 sq.deg. XMM-LSS survey l|Pierre et al.ll2007l ). 
Figure [l] presents the statistics of the exposure time dis- 
tributions. Figure [2] shows the distribution of hydrogen 
column density (Nh) on the corresponding lines of sight. 
Column densities values were obtained through th e Lei- 
den/Argentine/Bonn Nh maps (|Kalberla et alj |2005). Pro- 
cessed pointings show a median Nh of ~ 3.10 20 cm -2 and 
very few of them lie in regions above 10 21 cm -2 . We display 
on Fig. [3]the sky distribution of the processed observations. 
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Figure 3. Sky location of the 2774 pre-processed XMM observations (equatorial coordinates). Observations shorter than 10 ks (open 
circles) were not considered for the X-CLASS catalogue. The other pointings have been truncated to 10 ks (blue points) and, if possible, 
to 20 ks (gold points). Note that, by construction, a gold point (20 ks pointing) always has a 10 ks counterpart. The circle size is not 
representative of the XMM field of view (~ 30arcmin diameter). 



2.2 Data processing and cluster detection 

2.2.1 Processing steps 

The individual Observation Data Files (ODF) pertaining to 
each observation were retrieved via the ESA Archive Inter- 
operability system (XSAp]). Our processing is e ntirely based 
on the XMM-LSS pipeline l|Pacaud et al.ll2006l ) which main 
steps are summarized below. 

(i) Event lists are generated using XMM-SAS tasks 
emproc and epproc and filtered from proton and solar flares. 
This is achieved by creating the high energy events light 
curves (12 - 14keV for MOS and 10 - 12keV for PN) and 
flagging out periods of high event rates (rates greater than 
3 o~ above the mean observation count-rate). Although ad- 
equate for removing short periods of high flares, this pro- 
cedure may provide unsatisfactory results for observations 
having a high mean particle background. The overall qual- 
ity of each observation was then subsequently inspected by 
eye, and some observations discarded (Sect. [2~l3)l . 

(ii) Clusters detected with XMM exposure times of ~10- 
20 ks constitute the most relevant population for cosmolog- 

1 http://xmm.esac.esa.int/xsa/ 



ical studies l|Pierre et al.ll2007l . l201ll ): such exposure times 
are well above the XMM confusion limit. Further, the se- 
lection function of a survey consisting of homogeneous ex- 
posures is easier to handle. Using the good time intervals 
(GTI) resulting from the pre-processing, we thus build 10 ks 
and 20 ks "chunks" (now denominated "pointings" ) from the 
original exposures, starting at the beginning of the observa- 
tion. From each observation 0, 1 or 2 pointings are extracted, 
each pointing having exactly a 10 or 20 ks exposure time on 
the three detectors. The case of "0 pointing" occurs if one 
or more detector is insufficiently exposed, which corresponds 
to highly flared or problematic observations. Such observa- 
tions are discarded from the analysis, and the total number 
of pointings entering the source detection process is 2409 
(Table [TJ. We display in Fig. [3] the sky location of all point- 
ings having a duration of exactly 10 or 20 ks and that subse- 
quently undergo the source detection process. The archival 
GTI time used for cluster detection amounts thus to 24 Ms, 
over the 40 Ms GTI time available in total. In addition to 
these chunks, we construct a "full exposure" pointing con- 
taining the maximal GTI time from each observation, which 
will be used to obtain high signal-to-noise flux measurements 
for the detected clusters (Sect. |2~4")) . 

(iii) For each observation, three images are created in 
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Table 1. Number of XMM archival observations handled for the 
present study. First line, all available observations. Second line, 
retrieved observations. Third line, number of usable observations 
after event filtering. Last line, final number of retained observa- 
tions after discarding pathological cases (see. Fig. |S) 

XMM observations (May, 26th 2010) 7716 
Observations retrieved from the XSA 2774 

Pointings "chunks" 10 + 20 ks 2409 (1588 + 821) 

Pointings entering the cosmological analysis 1992 (1294 



the [0.5 - 2] keV band, for the three EPIC detectors, 
which are in turn co-added. The resulting image is sub- 
sequently filtered in the wave let space, assuming a Pois- 
son noise model (mr_f ilter, Starck. Murtagh. fc Biiaouil 



1 199Sl : IValtchanov. Pierre, fc Gastaudll200ll). and sources are 
extracted running SExtractor (|Bertin fc Arnoutsl Il996l ) 
over the filtered image. Only sources detected within the in- 
ner 13 arcmin radius FoV are considered in the subsequent 
analysis. 

(iv) Each detected source is further characterized by 
a maximum l ikelih ood profile fitting procedure (XAmin, 
iPacaud et all {2006)). Two source models are tested on each 
detection: (1) a point-like model using the analytical PSF 
from the XMM calibration database, with the position held 
at the SExtractor output value, thus allowing two param- 
eters to vary (the source count-rates on the MOS and PN 
detectors) ; and (2) a /3-model, convolved by the PSF, with 
five free parameters ((X,Y) position, core radius extent and 
MOS and PN count-rates). A uniform background is added, 
whose level is such that the total number of photons in the 
data equals that of t he modelled source plus background. 
We use the C-statisic jCashll 19791 ) for quantifying the likeli- 
hood of the fits and, finally, discriminating between the two 
types of sources. 

This procedure, when coupled to a representative set 
of simulations, readily enables the selection of sub-samples 
of X-ray extended sources with well characterised levels 
of completeness and contamination. It is thus perfectly 
suited to the analysis of large X-ray data sets such as the 
XMM archive. We note, however, that the XAmin pipeline 
was originally designed to detect and discriminate between 
point-like and exten ded sources in the XMM-LSS "empty" 
cosmological fields (|Pierre et al.l [2007); consequently, be- 
cause of the variety of astronomical objects present in the 
XMM archive (nearby galaxies, substructures in clusters, 
planets...) a subsequent human screening is necessary (see 
Sect. 12. 3|) in order to remove mis-interpreted detections. 



'point-like model' and is thus called the source Extent Like- 
lihood. These parameters can be easily related to the simu- 
lations intended to assess the survey selection function (see 
Sect. Because of their relevance, these values are listed 
in the final clust er catalogue (A pp.lAl). 

Following (IPacaud et al.l I2006T I. we denote by 
"CI" , sources characterised by EXTJLIKE > 33 and 
EXT > 5 arcsec : this corresponds to a sample of extended 
sources having a very low level of contamination by mis- 
classified point-sources. We further show in Sect. [3] that 
this parameter combination can be applied regardless of the 
pointing intrinsic properties. Fig. [4] displays the pipeline 
detections over one XMM archival pointing (ObsID: 
0403072201), containing three CI candidates and ~ 60 
point-like sources, most of them being AGNs. 



2.3 Catalogue construction 

2.3.1 Removal of duplicates 

Only the high-quality CI clusters are inserted into the final 
catalogue (App. \K\. Because of the multiple overlaps be- 
tween the archival observations, several sources are detected 
in more than one pointing. In particular, if an observation 
has been split in two pointings (10 and 20 ks), almost all 
sources detected on the 10 ks pointing are also found on the 
deeper one (26 over the 845 extended detections were not in 
this case, most of them being close to the detection limit and 
four being nearby, bright clusters with a very peaked profile, 
mis-classified as point-sources). Furthermore, because of the 
presence of the CCD gaps and/or of multiple maxima in the 
emission of widely extended sources, multiple detections of 
the same object occured. 

We first associated sources closer than 20 arcsec to each 
other. To decide which of the two sources has to be included 
in the final catalogue, we applied the following rules: If the 
two sources are on different pointings, the one lying on the 
deeper pointing prevails. If both detection lie on the same 
pointing or on different pointings having the same exposure 
time, the source with the lowest Extent Likelihood is dis- 
carded. Each association was controlled by eye to avoid the 
matching of close, unrelated sources (e.g. a cluster and a 
background AGN). Note that positional differences of a few 
tens of arcsec are possible for extended sources located at 
the edge of the FoV, mainly because part of the emission is 
missing and because of the distorted shape of the PSF. This 
procedure was thus repeated with larger correlation radii, 
until each catalogue entry was related to a unique source. 



2.2.2 Output parameters and source characteristics 

The procedures described above allow us to assign to each 
detected source a set of parameters characterising its prop- 
erties: position on the detector, off-axis distance, sky co- 
ordinates, count-rates in various energy bands, plus three 
numbers related to the chosen fitting algorithm; (a) the de- 
tection likelihood (DET_ML) gives the significance of the 
detection as compared to a pure background fluctuation; 
(b) the angular extent (EXT) is the apparent core radius 
of the best-fit /3-model; (c) the extent statistic (EXTXIKE) 
compares the significances of the 'extended model' and the 



2.3.2 Data screening and final selection 

All remaining entries underwent a detailed screening based 
on optical data. For each putative cluster, we retrieved im- 
ages from the Digitized Sky Survey (DSS) POSS-II on which 
we overlaid the X-ray contours. This step was mainly in- 
tended to remove extended sources, not relevant for our 
cluster catalogue: very large nearby clusters, halos of nearby 
galaxies, planets, unresolved double point-sources, and, pos- 
sibly, saturated point-sources. For this purpose, the DSS 
images are sufficient. During this process each source was 
assigned a quality flag by two astronomers independently; 
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Figure 4. Example of wavelet filtered image with detected 
sources (ObsID: 0403072201). 

the final decision was made by a moderator upon the evalua- 
tors' comments. False detections are classified as 'point-like', 
'double source', 'artefact' or 'nearby galaxy'. Among the 
1514 screened candidates, 234 X-ray detections were found 
to originate from nearby galaxies; 245 were classified as arte- 
facts, the majority of them being found in the X-ray emission 
of large, pointed galaxy clusters (Fig [5] bottom panels). An 
additional 'dubious' flag was assigned to sources for which 
the galaxy cluster nature is unclear: these mostly correspond 
to faint extended sources - at the CI limit - with an overall 
compact emission. As of Aug. 2010, the catalogue contains 
845 CI cluster candidates, 104 being classified as dubious 
(App.H) 

2.4 Count-rate measurements in multiple bands 

The XAmin pipeline is well suited to the count-rate mea- 
surement of moderately bright extended sources (up to a few 
thousands of photons). But in the case of sources occupying 
a large fraction of the detector or heavily contaminated by 
point-sources, hand-measurements are necessary in order to 
reach the accuracy required for an optimal application of the 
CR-HR method. We have thus developed a semi-interactive 
procedure to perform multi-band count-rate measurements 
and describe it below. 

2.4.1 Methodology 

We have developed software allowing for the masking of se- 
lected contaminating sources (mostly background or host 
AGNs), a careful account of the background levels and a 
possible redefinition of the source centre. Following paper I, 
we perform the measurements in the three different energy 
bands, well suited to the CR-HR analysis of the sample: 
[0.5-2], [1-2] and [0.5-1] keV. Input for the procedure are im- 
ages and exposure maps for the three detectors in the given 



bands. To correct for the masks, CCD gaps and detector bor- 
ders, the source to be measured is assumed to be spherically 
symmetric, and count-rates are integrated in concentric an- 
nuli. The initial X-ray center is the centroid determined by 
XAmin, but can be redefined by hand (e.g. if a CCD gap is 
masking part of the extended emission, and shifts the X-ray 
centroid). We define source count-rates as the mean num- 
ber of source photons collected by the telescope during one 
second, corrected from vignetting (i.e. equivalent to on-axis 
measurement) and detector cosmetics (CCD gaps, etc.). Us- 
ing the detector exposure maps, we derive the mean count- 
rate of the source in each annulus and compute uncertain- 
ties, assuming Poisson noise. A control annulus, chosen suffi- 
ciently far away from the source, but close enough to account 
for loca l variations, provides th e background estimate. Fol- 
lowing |Read fe Ponmanll2003t ). we model the background 
by a sum of a vignetted component (photon background) and 
an flat particle background, each of them being described 
by one parameter. The uncertainties on these parameters 
are derived assuming Poisson noise, and propagated to the 
uncertainties on the individual source count-rates. 

Measurements are performed on each detector, then 
summed up to provide a total count-rate. A count-rate 
growth curve is then computed (Fig. HU, as well as a surface 
brightness profile. For each source, the masking areas, the 
source position, as well as the number and width of the an- 
nuli and the background area, are set by hand in the [0.5-2] 
keV energy band. These settings are then stored and used for 
automatically measuring count-rates growth curves in the 
two other bands [0.5-1] and [1-2] keV. Final measurements 
are always performed using the complete pointing exposure 
(i.e. not only the 10 ks or 20 ks data), in order to maximise 
the signal-to-noise ratio. 

Such a procedure presents the advantage to be model- 
independent and does not require any spatial fitting. In turn, 
only aperture count-rates are available, up to a limiting 
radius at which the source emission vanishes in the back- 
ground. For each measurement, the software provides the ra- 
dius at which the integrated count-rate value shows a signal- 
to-noise ratio equals to 1 (i.e. compatible with background 
emission). The vertical dotted line displayed on Fig. [S]shows 
the position of this radius for a particular cluster measure- 
ment. In most cases, the integration radius set manually is 
very close to the automatic guess by the software, except in 
the few cases where the measurements in consecutive annuli 
are noisy (e.g. in presence of a mask or a CCD gap). 

2.4-2 Assessment of the method and aperture corrections 

We quantitatively evaluated the validity of our measure- 
ment procedure by means of simulated XMM observations 
of galaxy clusters. We used the simulation set described 
in Section [3] which provides a representative range of ob- 
serving conditions (in terms of background and exposure 
time) and of galaxy cluster sizes and fluxes. Table [5] sum- 
marizes the set of simulated observations. All simulated 
cluster profiles are azimuthaly symmetric /3-m odels with 
fixed P = 2/3 (|Cavaliere fc Fusco-Femiandl 19761 ). The total 
count-rate and core radius are taken among discrete values 
close to those expected in the survey. Each of the 18,000 
simulated observations is processed following the steps de- 
scribed previously. We then applied selection criteria identi- 



© 2002 RAS, MNRAS 000. 011251 



6 N. Clerc et al. 




ARC MINUTES ARC MINUTES 
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Figure 5. Illustration of the data screening classification (Sect. I2.3|l . From left to right, top to bottom : classification "0 < z < 0.3" 
(known cluster A2218 z=0.176), "z > 0.3" (indicative redshift, not used in the analysis), "Dubious", "Nearby galaxy" (NGC 4634). Both 
bottom panels show 'dubious' sources, likely being substructures in a close cluster or background clusters. Green circles indicate the 
position of CI candidates and green boxes are others X-ray detections; the cyan cross indicates the centroid of the X-ray emission; red 
crosses stand for the first-pass sources that were found to have a likelihood detection smaller than 15 by XAmin. 
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0.04 




Distance from center (orcsec) 

Figure 6. Example of a count-rate integration curve with associ- 
ated 1<t error bars (Cluster tag 2094 in L4SDB). The vertical line 
indicates a S/N of one. Apparent fluctuations at large cluster- 
centric distances are due to uncertainties in the background sub- 
traction. 

cal to those applied for the cosmological analysis, as outlined 
in Sect. 2] and the count-rate of the ~ 10, 000 selected clus- 
ters were automatically measured following the above proce- 
dure. The size of the background annulus and the boundary 
radius for the integration were set according to the cluster 
input extent. 

Fig. [7| presents the results obtained for both 10 ks and 
20 ks simulated pointings. Each panel shows the ratio of the 
measured count-rate over the true input value, as a function 
of measured core radii. The lack of statistics at large radii 
comes from the small number of extended sources detected 
in the simulations. The decreasing trend in this ratio as a 
function of input count-rate is explained by Eddington bias: 
at low fluxes, only clusters that pass the selection function 
criteria are measured. It artificially increases the mean value 
of the measured count-rate. Our method thus accounts for 
the statistical nature of the cluster sample. These results 
show the overall accuracy of the count-rate measurement 
procedure. They are in agreement with the fact that sources 
that are more extended are less well measured and this effect 
is more pronounced for higher background levels. 

To correct for the flux loss due to fixed aperture mea- 
surements, we fitted in each panel a linear relation of the 
form CRmos = tz-CRi llp ut + b (see Fig. [7]). We corrected indi- 
vidual cluster measurements by inverting this relation, thus 
providing an estimate of the true count-rate in the band of 
interest. 
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Table 2. Distribution of filter configurations for the set of XMM 
archival pointings entering the scientific analysis. Only 5% of 
those pointings have one or more detector observing with the 
Thick filter, which causes a ~22% diming at 1 keV. 

we form count-rate ratios obtained from a wide range of clus- 
ter spectra in the energy range of interest and fit a linear 
relation of the form: 

CR(X,Thin/Thin/Thin) _ / CR(F, Filterset) \ 

CR(X, Filterset) \ CR(X, Filterset) / ^ 

where X and Y are two different energy bands, the left-hand 
side representing the inverse of the filter attenuation in the 
considered energy band. We made use of Pimms and 
simulated XMM cluster count-rates on a grid of tempera- 
tures between 0.5 and 9keV and at different redshifts up to 
z — 1.5. Count-rates were computed in our three reference 
bands ([0.5-2], [0.5-1] and [1-2] keV). We did not consider 
clusters too cold and too dista nt since they are not retained 
by the CI selection function of lPacaud et all (|2006h . Fifteen 
filter combinations out of the 27 possible ones are found in 
the 1992 pointings used in our analysis (Tabled, more than 
half of them being in the ThinI-ThinI-ThinI configuration. 

We repeated the operation for typical galactic absorp- 
tion values, ranging from Nh = 10 20 to 2.10 21 cm -2 (see 
Fig. 01 ■ Table [3] shows an example of best- fitting values for 
the empirical correction. The corrections are significant only 
if the Thick filter is used, which is consistent with the fact 
that the ThinI and Medium attenuations are comparable 
for the chosen energy bands. The relative uncertainty of 
these conversions was found not to exceed a few percent. 



2-4-3 Filter combinations and Nh values 

According to the choice of the XMM guest observers, the 
MOS and PN observations are obstructed by one of the 
three EPIC filters, namely ThinI, Medium and Thick. For 
the purpose of applying the CR-HR method, cluster count- 
rates must be evaluated in a unique reference filter (THINI 
in our case). We thus need to apply some a posteriori cor- 
rection for using observations performed with other filters. 
Because filter transmission curves have different spectral de- 
pendences, these corrections are energy-dependent and we 
estimate them by means of empirical conversion relations: 



3 THE X-CLASS SURVEY SELECTION 
FUNCTION 

This section describes the steps leading to the construction 
of the final cluster sample used in the cosmological analy- 
sis. We present the list of retained XMM observations and 
subsequently describe the image simulations leading to the 
cluster selection function. We finally expose the method we 

2 http:/ /heasarc. nasa.gov/Tools/w3pimms. html 
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Figure 7. Aperture correction for measurements of cluster count-rates. Top: 10 ks simulated pointings, Bottom: 20 ks simulated pointings. 
In each panel only clusters having output core radii in the indicated range are included. Simulations were performed for various background 
levels as indicated by the colours. The data points have been horizontally shifted for clarity. Measured count-rates are obtained by the 
aperture photometry described in Sect. [2~4l Only simulated Cl + clusters (matching the same criteria as for the cosmological subsample) 
are included here. The black, dashed line shows the best linear fit used for correcting the individual count-rate measurements. 
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Table 3. Example of coefficients used for the empirical fil- 
ter corrections : bands X=[0.5-2]keV and Y=[0.5-l]keV and 
Njj = 4.10 20 cm -2 are considered here. The correction is modeled 
by a linear relation of the form x C orr = Cosmos + Ci .y mcs , where 
ic mcs and j/ mes are count-rates measured with the same particu- 
lar set of filters, respectively in bands X and Y and x COTT , the 
corrected count-rate as it would be measured in band X with the 
highest transmission set of filters (ThinI-ThinI-ThinI). These 
conversions were found to be accurate at the few percent level for 
clusters with temperatures ranging from 0.5 to 9 keV and redshifts 
out to 1.5, as long as their location in the temperature-redshift 
plane is covered by the CI selection (see text). 



have developed for evaluating the statistical bias due to the 
presence of numerous pointed clusters in the XMM archive. 



3.1 The cosmological subsample 

The cluster sample described in Sect. [2] contains 845 CI 
sources detected in a homogeneous way in the XMM archival 
data. In order to perform the cosmological analysis by means 
of the CR-HR method, we extracted a high signal-to-noise 
ratio subsample as follows. 

We first selected a more homogeneous data set, i.e. by 
excluding pointings (1) with a high background, (2) with one 
or more detectors not being in full-frame mode and (3) those 
centered on very nearby, luminous clusters (see examples on 
Fig- El)- For this purpose, we inspected the 2409 pointings by 
eye. In the end, the surveyed area used for the cosmological 
fits consists of 1992 pointings (Table [TJ) 

We defined a sub-class of galaxy clusters called Cl + 
by selecting all sources having an Extent Likelihood above 
40 and an extent greater than 5 arcsec. Only Cl + sources 
within 10 arcmin off-axis of their parent pointing and not 
flagged as dubious were considered. Defining CR as the [0.5- 
2] keV measured count-rate and HR as the ratio between the 
[1-2] keV and [0.5-] keV count-rates, we imposed clusters to 
have 0.009 < CR < 0.5 cts/s and 0.05 < HR < 2. Finally, 
the cosmological subsample consists of 347 clusters. 



3.2 Description of the simulations 

3.2.1 Principle 

We use an updated version o f InstSimulation 
jValtchanov. Pierre, fe Gastaudl l200ll ; iPacaud et al l 120061 ) 
to generate fake XMM observations, taking into account 
the main characteristics of the XMM EPIC instruments. 
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Figure 8. Three examples of pointings excluded from the cosmological analysis (images have been filtered by wavelets). Left panel: this 
pointing shows a high background on the PN detector which has not been optimally filtered by the pre filtering procedure due to long 
continuous periods of background flares (Obsld: 0039140101). Central panel: the PN detector was working in Large Window mode thus 
not collecting photons from the entire field of view (Obsld: 0083150401). Right panel /observation of Coma extending over a large part 
of the field of view thus preventing the detection of background sources (Obsld: 0300530301). 



Table 4. Number of detected clusters in the X— CLASS database. 
Clusters entering the cosmological analysis are Cl + sources not 
classified as 'dubious', within 10 arcmin off-axis and having 
0.009 < CR < 0.5cts/s (count-rate in [0.5-2] keV) and 0.05 < 
HR < 2 (ratio [1-2] keV / [0.5-1] keV). 

CI sources detected in the 2409 X-CLASS pointings 845 

C1+ sources : 745 

- classified as 'dubious', 74 

- within 10 arcmin off-axis, 630 

- entering the cosmological analysis 347 



In particular, an analytic vignetting model and a detector 
mask are superimposed to the simulated sources, as well 
as instrumental and photon backgrounds. The shape and 
off-axis dependency of the PSF is modeled by using the 
Medium model from the XMM calibration files. Apart from 
the peculiar observations of bright, saturated sources, or 
extremely extended sources, the simulation set captures the 
most important features impacting the detection and the 
characterization of the sources. We performed two sets of 
simulations: one with point-like sources only, and the other 
with clusters and point-like sources. 



3.2.2 Point-like sources only 

The first set of simulations without extended sources serves 
as a test for contamination and for parametrizing the back- 
ground level on the instruments. Point-like sources are 
distributed across the fi eld of view, followin g a sampled 
logA^-logS 1 taken from (|Moretti et alj 120031 ) in the [0.5- 
2 keV] band. Conversion from flux to count-rate is performed 
assuming a constant MOS to PN count-rate ratio, regardless 
of the source spectral distribution. The flux lower bound is 
chosen accordingly to the exposure time so as to give ~ 2 
photons on-axis (i.e. below the XMM detection limit). Non- 
resolved AGN photon ba ckground is added following val- 
ues from (Rea d fc Ponmanll2003h . then corrected from the 
estimated fraction of AGNs resolved by the pipeline. This 
background component is vignetted, thus showing a strong 



off-axis dependence. We finally add the non-vignetted par- 
ticle background c ompo nent using the standard values from 
I Read Sz Ponmanl |2003h multiplied by an arbitrary factor 
6 = 0.25, 0.5, 1, 2 or 4 so as to investigate the impact of 
pointing-to-pointing background variations on the detection 
efficiency. For each background value, 540 pointings are sim- 
ulated both at 10 and 20 ks exposure times. 

Each of these pointings is processed by the XAmin 
pipeline described in Sect. [2j exactly in the same way as 
real observations. In particular, an extended model fit is 
performed over each detected point-source to evaluate the 
contamination of the cluster sample. 



3.2.3 Extended source simulations 

Similar simulations are performed by adding clusters as- 
sumed to be spherically symmetric sources and defined by 
a /3 = 2/3 profile; apparent core radii r c range from 10 to 
100 arcsec and total count-rates from 2.5. 10 -3 to 0.1 cts/s. 
From 4 to 8 simulated clusters are injected in each point- 
ing, depending on their angular size, and we avoided source 
overlap by defining exclusion sectors in the XMM field of 
view. In total, some 87000 extended sources were simulated 
over more than 18000 pointings (Table [SJ). In order to repro- 
duce the effects of point source contamination, a population 
of point-like sources was added as described in the previous 
section. 

Each pointing is processed by the pipeline and positions 
of the input extended sources are correlated with the out- 
put positions within a 37 arcsec radius. In case of multiple 
matches, the detected source with the highest Extent Like- 
lihood is chosen as the best matching counterpart and all 
others are discarded. Fig. (top panel) shows three exam- 
ples of simulated pointings in different observing conditions, 
along with the corresponding pipeline results. 
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Table 5. Summary of the extended-source simulations in XMM 
images. Last column indicates the number of simulated clusters 
out to an off-axis of 10 arcmin. Simulations were performed for 
exposure times of 10 and 20 ks and five background levels (b=0.25, 
0.5, 1, 2 and 4). The total number of simulated and processed 
pointings is 18140 and the total number of clusters amounts to 
~ 87000. 



Input count-rate 


Input core radius 




(10~ 2 cts/s) 


10" 


20" 


50" 


100" 


Total 


0.25 


450 


450 


750 


240 


1890 


0.5 


450 


450 


750 


240 


1890 


1.0 


450 


450 


750 


240 


1890 


2.0 


450 


450 


750 


240 


1890 


5.0 


100 


100 


150 


240 


590 


10.0 


100 


100 


150 


240 


590 



3.3 Analysis of the simulations - selection criteria 

3.3.1 Contamination by spurious and point-like sources 

Following iPacaud et al.l (|2006l ) we report for each detected 
source its location in a two-parameter space Extent/Extent- 
Likelihood, as shown in Fig. [9] middle panels. In this figure, 
green symbols represent point-like sources, magenta sym- 
bols are for extended sources and red points stand for spuri- 
ous detections, i.e. detections in the point-only simulations 
which are not associated to input sources within the 6" cor- 
relation radius. Fig. [9] shows the good stability of the Cl + 
criterion across the range of exposure times and background 
levels, in terms of contamination by point-like sources and 
spurious detections. From our simulations, we expect the 
number of contaminating sources (i.e. point-sources inter- 
preted as extended sources) not to exceed one in every 300 
pointings for normal observing conditions (10 ks, low back- 
ground) . 



3.3.2 Efficiency of the extended sources detection 

We then derive the Cl + detection efficiency by taking all 
sources in the 0-10 arcmin off-axis range. Bottom panels 
of Fig. [9] display the probability of detecting a Cl + clus- 
ter as a function of its total input count-rate and its input 
core radius, as derived from our simulations. From these re- 
sults, it clearly appears that the selection is not flux-limited, 
but rather surface-brightness limited. These curves also re- 
veal the expected increase in efficiency from 10 to 20 ks and 
for lower background levels. The sharp decrease observed in 
the 20 ks selection for high count-rates (~0. 06-0.1 cts/s) and 
small core radii (< 20 arcsec) indicates that these sources are 
identified as point-like sources by the pipeline; such objects, 
however, are unlikely to appear in real observations. We de- 
rived similar probability functions for the ten simulated con- 
figurations (10 and 20 ks exposures and five background lev- 
els). 

Then, for any given pointing in the survey, we estimate 
its background parameter (b). To this purpose, we measure 
local background estimates at several locations on the de- 
tectors and compare them to the values found in the sim- 
ulated observations for which b is known. We finally derive 
the selection function of each XMM observation entering the 
cosmological analysis. 



3.4 The survey geometrical area 

For the cosmological analysis, we consider only sources 
within a 10 arcmin radius around each pointing center, thus 
removing regions of the detectors where the point spread 
function has an elongated shape and the vignetting factor 
is greater than 50%. Because of the multiple overlaps be- 
tween pointings, we estimate the net area by means of a 
Monte-Carlo integration. If two or more pointings of the 
same exposure time overlap, the intersecting area is equally 
distributed between those pointings. If one 20 ks pointing 
and one 10 ks pointing overlap, the intersecting area is fully 
attributed to the deeper pointing. This set of rules is thus 
compatible with the procedure applied for removing dupli- 
cate detections in the catalogue (Sect. 01. The net total area 
covered by the 1992 pointings is 90.3 deg 2 . 



3.5 Correcting for the target bias in the XMM 
archive 

In all cluster serendipitous surveys to date, it has always 
been implicitly assumed that discarding the central target of 
the considered pointings (along with subtracting the corre- 
sponding survey area) would not statistically affect the stud- 
ied cluster population. It is not the purpose of the present 
paper to verify these past assumptions, but in the case of 
serendipitous surveys based on the XMM archive, this hy- 
pothesis is questionable. Among the 7716 archival observa- 
tions available by May 2010, 1008 of them pertain to pointed 
observations of galaxy clusters. Out of the 347 clusters se- 
lected for the present cosmological analysis, 92 of them are 
central targets (within 3 arcmin offaxis). One cannot simply 
discard them or include them (or ignore the complete point- 
ing) in the statistical analysis, because the process of target 
selection from guest observers is highly subjective as well as 
motivated by practical constraints. This is particularly true 
for distant clusters, as only the brightest ones could be ob- 
served. 

In this section, we describe the method that we have devel- 
oped to account for the target selection bias and correct for 
its impact in the cosmological analysis. This bias is hard to 
model from first principles as it depends on the history of 
the XMM observing programs. Basically, we make use of the 
fact that pointed galaxy clusters are preferentially located 
at the center of the XMM field of view and split our sample 
in two subsamples to apply a joint correction. 



3.5.1 Off-axis source distribution 

Fig. [10] shows the distribution of detected clusters as a func- 
tion of off-axis distance. The excess of sources in the 0- 
5 arcmin range is conspicuous, as well as the stabilization at 
higher radii. There is a factor of 2 between the density of 
sources in the 'outer' regions of the detectors and the 'in- 
ner' regions. Part of this difference stems from the better 
sensitivity in the central part of the XMM FoV but a also 
from a number of pointed clusters. We label by 'in' and 'out' 
the corresponding two subsamples and make the hypothesis 
that all pointed clusters are found in the 'in' sample. 
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Figure 9. The X-CLASS selection function. We present the results of simulations performed for three pointing configurations (different 
exposure times and background levels) among the 10 simulated configurations. Top: Example of simulated XMM observations. Blue 
circles show the locations of simulated 20arcsec core-radius clusters and magenta circles denote detections classified as Cl + by the 
pipeline. Green boxes are for the remaining, unclassified sources (including point sources). The radius of the large blue circle is 13arcmin. 
Middle: Distribution of detected sources in the Extent-Likelihood/Extent plane. Recovered clusters are in magenta, point-like sources 
in green and spurious detections in red. Vertical and horizontal lines delimit the C1+ classification (EXT_ML > 40 and EXT>5); the 
discreetness of the Extent distribution reflects the input core radius values (Table |5j Bottom: Corresponding detection probability for 
the CI"*" sources as a function of the input [0.5-2] keV count-rate and the input core-radius. [2. 5. 10 -3 , 5. 10 -3 , 0.01, 0.02, 0.05, 0.1] 



3.5.2 Bias model 

We display in upper panel of Fig. QT] the count-rate dis- 
tribution of sources in both subsamples per unit sky area. 
According to Fig. 1101 there is a factor ~2 between both 
distributions, but the the excess is not uniformly spread 
over the range of count-rates. For instance the factor 3 ex- 
cess around CR~ 0.15 cts/s corresponds to clusters with 
fluxes ~ 1 — 2.10 -13 ergs/cm 2 /s typical of those found 
in ROSAT serendipitous surveys: 160d (|Vikhlinin et al.l 



1998b), 400d dBurenin et all |2007h. WARPS (I Jones et al l 
19981 ) and SHARC (lAdami et a.1.1 l200ol : iRomer et alJl200Ch 
- see iPiffaretti et all (|201ll ) for a thorough compilation of 
ROSAT cluster catalogues. Clusters from the ROSAT All- 
Sky Survey are more than ten times brighter on average and 
are thus excluded from our cosmological sample (limited to 
CR < 0.5 cts/s). 

Because of the finite number of clusters of given flux 
across the entire sky, the sample dubbed 'out' does not ex- 
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Figure 10. Surface density of clusters as a function of off-axis 
distance, for the selected subsample used in the cosmological anal- 
ysis (black curve). The red and blue points are the average sur- 
face densities in the inner [0-5] arcmin and the outer [5-10] arcmin 
regions respectively. The factor ~ 2 between the two values is 
mainly caused by the presence of pointed clusters in the archival 
data. 



actly reflect the cluster population as it lacks all clusters 
being pointed. As both subsamples derive from the same 
parent distribution, we use a single parameter for the inner 
excess and the outer dearth of clusters. We detail in App.lBl 
our procedure to infer its value, taking into account the ef- 
fective areas of both the inner and outer part of the XMM 
field of views. All count-rate bins j are treated separately 
and we compute a bias factor Fj (Fig. I12[) whose value rep- 
resents the ratio between the observed number of clusters in 
the considered bin and the actual expected number of clus- 
ters if no object were pointed. By definition, Fj is always 
greater or equal to one. 



4 THE COSMOLOGICAL ANALYSIS 

This section presents the analysis performed with the sub- 
sample of 347 Cl + clusters selected over our effective area of 
90 deg 2 from the XMM archival data. We first show the re- 
sulting CR-HR diagram, which is the sole observable quan- 
tity used in the cosmological analysis. We then describe its 
modeling from first principles, taking into account a cos- 
mological model, X-ray cluster scaling-laws and the various 
selection effects affecting the sample. We finally show the 
results obtained by a Monte-Carlo Markov Chain sampling. 



4.1 Sample CR— HR distribution 

Following Sect.O count-rates for all 347 clusters entering the 
analysis have been measured in three energy bands: [0.5- 
2], [0.5-1] and [1-2] keV. These values have been corrected 
from flux loss due to the finite aperture measurement and, 
if necessary, converted in the ThinI filter set. We compute 
the hardness ratio of each cluster by dividing the [1-2] keV 
count-rate by the [0.5-1] keV measurement and report its 
value in a CR-HR plane, where CR stands for the wide 
band measurement. Fig. [13] shows the distribution of clus- 
ters in this diagram, along with associated error bars. The 
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Figure 11. Top panel: count-rate distribution of the sources 
found within 5 armin off-axis (red) and between 5 and 10 ar- 
cmin off-axis (blue) for the cosmological subsample containing 
347 clusters (all count-rates are rescaled to their on-axis values). 
Bottom panel: the plain curve shows the ratio between the two 
histograms shown in the top panel. The dashed line is the ratio 
that one would expect from the sensitivity gradient only. 



distribution is spread over the range of count-rates peaking 
around 0.3-0.4 cts/s. 



4.2 Modelling the CR-HR distribution of sources 

We describe the main steps of the computation of the CR- 
HR distribution of clusters, starting from the halo mass func- 
tion and using the survey selection function. These steps are 
more thoroughly detailed in paper I. 



4-2.1 Halo mass and redshift distribution 

We assume a ACDM cosmological model and a flat Universe 
(f2 m + S7a = 1) with no-evolving dark energy (w — — 1). 
Starting from the primordia l, scale-invariant po wer spec- 
trum with slope n s = 0.961 (iDunklev et al.l [2009l l we make 
use of the lEisenstein fc HvJ ( 19981 ) fitting formula for the 
transfer function to obtain the linear power-spectrum. We 
use the iTinker et all l|200Sl ) fit to the mass function and 
obtain the comoving density of haloes per mass interval 
dM2oob about M2oob at redshift z, where M2oob is the mass 
within a radius i^oob inside which the mean density is 200 
times the mean density of the Universe at that redshift. 
We convert this distribution into the sky-projected density 
of clusters per redshift slice. Only f2 m and as (normaliza- 
tion of the power spectrum at a scale R = 8ft" 1 Mpc) are 
let free in the analysis . All other paramet ers are held at 
their WMAP-5 value (|Dunklev et all I2009T ), in particular 
H = 72.4km/s/Mpc. 
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Figure 12. Bias factor Fj for cluster counts (see App [B} due 
to the presence of pointed clusters in XMM archival data. Fj is 
computed in 10 count-rate bins and evaluated by comparing the 
density of clusters in the inner [0-5] arcmin to the [5-10] arcmin 
density, taking into account the higher sensitivity around the 
EPIC optical center. The bias value is the ratio between the ob- 
served number of sources in a given bin and the number expected 
for a sample without pointed clusters. We use it as an empiri- 
cal correction for the final sample of 347 clusters selected for the 
cosmological analysis. 



0.01 0.10 
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To account for the intrinsic scatter in those relations, we as- 
sume two parameters <r inT | A/ and o~i nL \ T constant over the 
entire redshift, mass and temperature ranges considered in 
our analysis. T hroughout this work we will use the M2 00C- 
T relation from lArnaud. Pointecouteau. fc Pratt! |2005l ) de- 
rived for their hot cluster sample (see Tabled]). 

As discussed in paper I, a reasonable choice 
for the emissivity profile of clusters is a f3- model 
(|Cavaliere fc Fusco-Femianol |l976) with /3 = 2/3 and a 
core radius r c scaling with Tlsooc, parametrized by a; c ,o = 
rc/Rsooc at all redshifts and masses. This parameter is crit- 
ical as it enters in the selection function describing the ob- 
served population of clusters. 



4-2.3 Instrumental model and measurement errors 

The next steps consist in converting z, T and Lx into ob- 
servable quantities and fold the cluster distribution into the 
selection function in order to obtain a CR-HR diagram for 
a given set of model parameters. Ideally, one should com- 
pute this distribution separately on each pointing and sum 
up their independent contributions to derive the complete 
CR-HR diagram. In order to avoid excessively large compu- 
tational times we decided to group pointings by similar Nh 
and background values. 

We use APEC spectral models with a metallicity of 
0.3 Zq along with EPIC XMM response matrices to obtain 
count-rates in the three bands of interest. They represent 
the expected number of source events per second impact- 
ing the telescope cameras mounted with ThinI filters. The 
[0.5-2] keV count-rate as well as the cluster apparent core- 
radii are then derived. Finally, the sample selection function 
is used to compute the expected distribution of clusters. 
Measurement errors are included by convolving the CR-HR 
distribution with an error model using the statistical uncer- 
tainties from the count-rate measurements. 

To fully model the cluster population we add a supple- 
mentary step to the methodology presented in paper I: the 
CR-HR model distribution is multiplied by the bias value 
F(CR) as computed in Sect. [3] (see Fig. I12[) to account for 
the excess of sources due to pointed observations. 



Figure 13. CR-HR diagramme for the 347 clusters pertaining 
to the cosmological subsamplc. Not all error bars are displayed in 
order to ease visualization. The dashed box delimitates the region 
used for the cosmological fit. 



4-2.2 Cluster emissivity and extent 

We assign to each cluster three quantities describing its X- 
ray emission: its redshift z, its plasma temperature T (X- 
ray spectral temperature) and its bolometric luminosity Lx 
integrated over its complete extent. We consider a mean 
metallicity of 0.3 Zq. The conversion from cluster mass to 
temperature and luminosity is achieved thanks to scaling 
relations of the form: 



A-f200c 



10 



(ikiv) E (*r\i + *r^ (2) 
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4.3 Likelihood and MCMC sampling method 

Given a set of parameters, the expressi on for the likelihood L 
is expressed following e.g. ICashl (|l979h by dividing the CR- 
HR two-dimensional space in narrow bins such that each bin 
contains at most one cluster. Assuming Poisson statistics in 
each bin, we can write: 

dn 



ln£ = 



dCRdHR 



CR„ 



(CRi,HRi) 



dn 
dCRdHR 



dCRdHR 



(4) 



where we have neglected the constant term including the 
size of the bins as we will consider likelihood ratios only. 
In the equation above, the sum runs over the 347 selected 
clusters. The integral is computed for 0.09 CR ^ 0.5 
and 0.05 ^ HR ^ 2 and simply represents the expected 
number of clusters within this CR-HR region. The modeled 
number density contains the effect of measurement errors 
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and pointed cluster bias and includes the selection function 
and its variations over the surveyed area, thus is as close 
as possible to the real CR-HR distribution of sources. We 
summarize in Table [6] the choice of parameters and priors 
we made for the different cases studied. 

Confidence intervals and mean values for the parame- 
ters being studied are obtained via the Bayesian formalism 
and are computed from the posterior distribution of parame- 
ters given the CR-HR diagram and the underlying model. In 
this study, we limit ourselves to a maximum of 5 free param- 
eters, for which a Monte-Carlo Markov Chain (MCMC) like- 
lihood exploration becomes competitive with a grid-based 
computation. We use a custom Metropolis-Hastings sam- 
pler that produces MCMC chains. After a so-called burn-in 
period, the chain reaches a stationary state representative 
of the actual posterior parameters distribution. The 'jump' 
function is taken as a multivariate gaussian distribution with 
covariance matrix C M „ computed with (see paper I for a more 
thorough description and references therein): 

where Oi stands for the binned density dn/dCR/dHR and 6 
is the set of varying parameters. The PDF 'jump' function 
from the current parameter set n to the new one thus 
writes : 

p{0\0 n ) ocexp[- i T (6l-6» n )F^(6»-6> rl )] (6) 

For each chain we extract the mean value of the sam- 
pled parameters and compute the associated highest den- 
sity intervals. Such an interval contains (1 — a) 100% of the 
posterior probability and ensures that the posterior den- 
sity within the interval is always greater than outside. We 
choose a — 0.32 which, in the case of a normal distribu- 
tion, corresponds to the 1-a boundaries of the distribution. 
The MCMC analysis is performed using the R-package BOA 
l|Smith|[2007l ). Unless otherwise stated, best-fit results are 
quoted using the mean of the posterior distributions (and 
not the maximum likelihood estimate). 



5 RESULTS 

Our methodology, as presented in paper I, assumes that 
local scaling laws are known, we thus focus on their evo- 
lution and on the cosmological parameter determination. 
Different expressions for the local scaling laws are found 
in the current lite rature. This is particularly true for the 
L-T r elation (e.g. lArnaud &: Evrard| [T999; Branchcs i et al.l 
120071 : iMaughan et all 120121 ), likely beca use of the differ- 
ent populations b eing taken into account l|Pratt et al.ll2009l : 
iMittal et al.ll201lh . but also because of different selection ef- 
fects. However, our cluster sample is not quite large enough 
to allow for a simultaneous fit of the local scaling laws, of 
their evolution and of cosmology. We thus proceed with a 
step by step approach to select a local L-T t hat matches 
well o ur sample, focusing on two relations from lPratt et al.1 
(2009). We use these relations since they are well suited to 
the ma ss range of our sample, contrary to, e.g. iMantz et al.l 
l|2010bh scaling relations which have been derived for much 
more massive clusters. We first set the cosmology to the 



WMAP-5 values and select a 'best' relation by comparison 
to our data and to published logTV-logS from the literature. 
We then release cosmological parameters to perform an en- 
larged fit of our data and finally show that the selected L-T 
still adequately describes our clusters. In all cases, we always 
assume a flat ACDM Universe. 

5.1 Fixed cosmology, fixed local scaling laws 

5.1.1 Constraints from the CR-HR distribution 

We first fit 7 z ,mt and 7z,lt, the parameters governing the 
non-self-similar behaviour of the M-T and L-T relations, as 
well as x c ,o, the cluster size. Other parameters are set ac- 
cording to Table [6j in particular Q m and as which are held 
at their WMAP-5 value and we assume that scaling rela- 
tions are perfectly known. For the temperature t o luminosity 
conve rsion we consider the numerical values of IPratt et al.l 
(2009) using two of their L\ — T\ relation without core exci- 
sion. One of them has been computed for all clusters present 
in their sample ('ALL') while the other ('NCC') excludes all 
clusters showing high central gas density and thus hosting a 
cool core. As the intrinsic scatter in these relations depends 
on the population of clusters in the sample, we use two test 
values, oi n L\T = 0.3 and 0.7; these valu es correspond to th e 
lowest and the highest scatter found in lPratt et al.l ([2009) . 
We emphasise here that we allow us to somewhat generalise 
their results since, formally, scatter values of 0.3 and 0.7 are 
associated to the 'NCC and 'ALL' samples respectively. 

The fit results with associated uncertainties are quoted 
in Table [7] along with an indication of the relative likelihood 
for each test case. The Cash statistics C = — 21n£ is evalu- 
ated at the mean of the posterior distribution. As changing 
from one L-T relation to the other involves 3 parameters 
in our model («lt, Clt, 7z,lt), the difference between C 
and the (unknown) minimal value C m in obtained when all 
6 param eters are let free behaves as a \ 2 with 3 degrees of 
freedom (|Cashlll979h . This allows us to put a lower boundary 
to the x 2 of au four fits, e.g.: 

Xall,o.7 — Call, 0.7 — C m m 

> Call, 0.7 — Cncc,o.7 = 8.4 

Because the probability for x\ t° De greater than 8 is 
~ 5%, it indicates that our data prefers the 'NCC scaling 
relation with a large intrinsic scatter to the three other scal- 
ing laws. 

5.1.2 Comparison with published logN-logS distributions 

We cross-check the outcome of our 3-parameter fit by com- 
paring the resulting cluster logA-logS 1 to that from other 
surveys. We show on Figure [14] the logiV-logS computed 
with two local L-T ('NCC and 'ALL', a lnh \ T = 0.7). We 
display the result with and without the (1+z) 7 evolution fac- 
tor in the scaling laws. In the former case, we used our best- 
fit parameters 7 z ,mt and 7z,lt from Table [7] We note that 
the value of x Ct o has no impact on the predicted logA'-logS 
since such a distribution is meant to be flux l imited . The 
REFLEX data correspond to iBohringer et al.1 (120011 ) best- 
fit power-law converted into a [0.5-2] keV logiV-logS' using a 
constant factor calculated for a AP EC plasma at z = with 
T = 5 keV. RDCS values are from iRosati et ail (|l998T ) and 
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Table 6. List of parameters used in this work. The cosmological parameters are from WMAP-5 JPunklev et alJl200Sh . Numbers in 
brackets indicate the uniform priors applied in the MCMC fitting procedure. When fixed, f2 m and as are held at their W MAP-5 values, 
namely O m = 0.249 and as = 0.787. 'ALL' and 'NCC refer to the corresponding L— T relations from lPratt et all j2009f ). 



Parameter 


Fixed value or [prior range] 


Description 




[0.09 - f] 




n A 


1 - fin, 


(Flat Universe) 


J " b 


0.043 




CT8 


[0.05 - 2] 




n s 


0.961 




h 


0.72 






1.49 


Af — T power-law index 


Cmt 


0.46 


M — T logarithmic normalization 


7z,MT 


[-4, 4] 


M — T evolution index 


a lnT\M 


0.1 


M — T constant logarithmic dispersion 


a LT 


'ALL': 2.7, 'NCC':2.9 


L — T power-law index 


Cut 


'ALL': 0.52, 'NCC: 0.40 


L — T logarithmic normalization 


7z,LT 


[-5,3] 


L — T evolution index 


u lnL\T 


0.3 or 0.7 


L — T constant logarithmic dispersion 


Xc,0 


[0 - 0.9] 


/3-model core radius scaling wrt. iJsooc 



Table 7. Best-fit values for the evolutionary parameters 7 z ,mt, 7z,lt and the geometrical scaling factor x c ,o for a fixed cosmology. 
Quoted results are the mean and 68% confidence intervals obtained by fitting the X-CLASS CR-HR distribution, while cosmological 
parameters are held fixed at their WMAP-5 value and only those 3 parameters are varied. The C = — 2ln C m ax values are computed at 
the location of the best-fit parameters and differences are quot ed relative to that obtained for 'NCC with large scatter (2nd column). 
'ALL' and 'NCC refer to the corresponding L— T relations from lPratt et all ll2009h 



Local L-T: 



NCC 



ALL 



a lnL\T'- 


0.3 


0.7 


0.3 


0.7 


7z,MT 


0.60 ±0.15 


0.32 ±0.13 


0.17 ±0.18 


-0 1S+ 016 


7z,LT 


-1.23 ±0.41 


-1 30+ 054 


9 9 <-+0.61 
z - zd -0.48 


-2 06 +0 - 56 




0.17 ±0.02 


0.26 ±0.03 


0.27 ±0.02 


0.39 ± 0.04 


C — CnCC,0.7 


7.9 


(0) 


9.1 


8.4 



160d a nd 4 00d correspond respe ctively to IVikhlinin et al.l 
l|l998bh and lBurenin et alJ (|2007r ). 

From this figure it turns out that the 'ALL' relation 
predicts too many clusters compared to observations from 
other authors, particularly for the brightest, most massive, 
nearby clusters. Conversely, the 'NCC scaling relation (also 
assuming o lnL \ T = 0.7) is less discordant and our proposed 
evolution nicely fits the low-flux logA r -log5 from the RDCS 
and the 400d surveys. 



5.2 Free fi m and as, fixed local scaling law 

We now relax S7 m and as while fitting the CR-HR diagram, 
in addition to 7z,mt, 7z,lt and x c ,o, thus allowing 5 param- 
eters to vary in total. We consider the 'NCC relation with 
a inL\T — 0.7 as our reference L-T relation and keep it fixed 
in the analysis. Figure [15] shows the resulting posterior dis- 
tribution obtained from the MCMC chains. The mean value 



and associated 1-a error bars for each parameter are : 

o - n 94+ - 04 

i'm — u-^t-o .09, 

as = O.88I0.13, 

n oq+0.45 
7z,MT = U.«.3_o.56, 

7z,LT = -l.d_ .7> 

x Cj0 = 0.24 ± 0.04. 

We note that the values for the three last parameters are 
consistent with the previous results from the 3-parameter 
fit. Figure [16] illustrates the good agreement between the 
observed CR-HR distribution and the best-fit model. Using 
the best-fit model we predict a total amount of 369 clusters 
in the region where the fit is performed (red dashed box 
on Fig I16|) which is comparable to the 347 clusters actually 
present in the analysis. We note the presence of a few outliers 
that will be discussed in the next section. 

We show on Figure [17] the cluster logA r -log5' predicted 
by this set of best-fit parameters, along with the 1-cr lower- 
and upper-boundaries obtained by propagating the posterior 
covariance matrix from the MCMC analysis. The logTV-logS 
agrees well with the reference logTV-logS curves, especially in 



the range 10 



2.10 ergs/s/cm . This interval roughly 



corresponds to the count-rate region probed by our CR- 
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Figure 14. The cluster logTV-logS 1 . The blue and red curves correspond to our predictions assuming the local 'NCC and 'ALL' L— T 
relations respectively, both taken with <r lnL \ T = 0.7 and are computed for the WMAP-5 cosmology. Left panel: Self-similar evolution is 
assumed for the scaling laws (7 Z kt = 7z lt =0)- Right panel: Evolution as inferred from our 3- parameter best fit on the X— CLASS 
CR-HR distribution (7 Z mt an d 7z LTi Table [TJl. Data points correspond to observations from various surveys (see text). The vertical 
lines indicate the approximate flux range of our sample. 



«j B%aa S a , „ J '□ % " 

«Sn°¥ ^ if S " =» " 



1 



0.01 0.10 

[0.5-2] keV count-rote (CR, cts/s) 




1x10"* 
10 



10-'" 10~ 13 10 - ' 2 

S (ergs/cm 2 /s) in [0.5-2] keV 



10-" 



10-'° 



Figure 16. Best fit model (green dashed contours) overlaid on 
the CR-HR diagram data points. Five parameters are fitted as- 
suming the local L-T relation from 'NCC and c ini | T = 0.7. 
Contour labels indicate the predicted number of clusters enclosed 
by the corresponding contours. The red box shows the region in 
which the fit is performed. Measurement errors and bias due to 
pointed clusters have been included in the computation of the 
green contours. 



Figure 17. The cluster logW-logS. The plain red curve shows 
the modeled distribution computed assuming the local L— T rela- 
tion of 'NCC with cr ;ni | T = 0.66. Cosmological parameters as 
well as the non-similar evolution of scaling laws have been ad- 
justed so as to match the CR-HR distribution of sources in the 
X-CLASS sample (best-fit parameters quoted in Sect. [52}. The 
red dotted curves are computed by propagating the uncertainties 
on these parameters in the log./V-logS computation. The vertical 
lines indicate the approximate flux range of our sample. 



6 DISCUSSION 



HR diagram (assuming an average flux conversion factor of 
9.1CP 13 ergs/s/cm -2 per cts/s). The high- and low-flux ends 
of this curve are not probed by our data points but rather 
rely on the validity of the assumed model, and in particu- 
lar on the fact that scaling laws behave as Equs. [2] and [3] 
at all redshifts. Such an extrapolation disagrees with RE- 
FLEX data points, as is the case in Fig. [14] where cosmo- 
logical parameters are held at their WMAP value. Changes 
in the shape of scaling laws as a function of cluster prop- 
erties and/or a disagreement between the different selection 
functions of various surveys may explain this discrepancy. 



6.1 Cosmological parameters 

Our analysis indicates values for f2 m and erg of 0.24 and 0.88 
respectively with a ~ 25% and ~ 15% accuracy. These con- 
straints are co mpatible with the most recent measurem ents 
from the CMB jDunklev et al.ll2009l;lLarson et al.ll201ll) an d 
BAO observations (|Percival et all |2010| ; iBlake et all 1201 ll ). 
They are also in agreem ent with the most rec e nt studies of 
X-ray selected clusters ( Vikhlinin et al.l 120091 ; Irlenrv et al.l 
120091 ; iMantz et al.ll2010ah. Sunyaev-Zeldovich sel ected clus- 
ters (|Vanderlinde et all 120101: ISehgal et al.|[201ll ) and opti- 
cally selected clusters l|Rozo et alj 20ld ). 

Our results have been obtained for fixed local scaling 
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Figure 15. Posterior distribution for the five parameters fitted using the X— CLASS CR-HR distribution. The local L-T relation is 
'NCC or 'ALL' with a logarithmic dispersion <r ini | T = 0.7. Diagonal panels represent the one-dimensional marginal distributions 
for each parameter (normalized to unit area), sub-diagonal panels show the two-dimensional contours enclosing 68% and 95% of the 
marginalized posterior distribution. 



relati ons, in part i cular the 'NCC L-T relation was taken 
from IPratt et al.l l|2009h . with a constant logarithmic scat- 
ter of 0.7. We note that their relation has been derived for an 
assumed ACDM cosmology (fi m = 0.3,fiA = 0.7,h = 0.70) 
and strictly speaking it should be converted for each tested 
cosmology. However, we checked that the correponding cor- 
rection from f2 m = 0.3 to fi m = 0.24 on the L-T normal- 
isation amounts to less than 5% and thus neglected this 
correction in the analysis. 

We checked that assuming the L-T relation from 



lArnaud fc Evrardl jl99$ ) (converted in the WMAP-5 cos- 
mology) leads to compatible constraints on f2 m and as- In - 
terestingly, the 'ALL' (a [nL]T = 0.7) from lPratt et all (|2009f) 
predicts fi m = 0.15 ± 0.04 and erg = 0.96 ± 0.08, in disagre- 
ment with the other works cited above. Fig compares the 
confidence contours obtained for both scaling relations and 
illustrates the marginal agreement of t he posterior distri- 
butions. In all cases ('NCC, 'ALL' and lArnaud fc Evrardl 
1999) we tested two different values for the intrinsic scatter 
&inL\T (0.3 and 0.7) and noticed little change in the (fi 
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Figure 18. Predicted redshift distribution for the "01"'"" clusters 
in the X-CLASS cosmological subsample (red solid histogram). 
The redshift distribution of clusters in our sample with known 
redshifts from the literature is indicated by the black dashed line. 



constr aints. This result is compatible with ISahlen et al.l 
(2009) who find very little degeneracy between <7 ini | T and 
cosmological parameters. 

From the best set of parameters (Sect 15. 2[) . we infer 
the redshift distribution of our cosmological cluster sam- 
ple (Fig. I18|) for the pure Cl + selection (with no bias from 
pointed clusters). It shows (a posteriori) that the distribu- 
tion peaks around z me d ~ 0.3 — 0.4 with only 8.4 Cl + clus- 
ters beyond z = 1. The redshift histogram of clusters with 
known redshifts from the literature clearly illustrates the 
complexity of this bias in our sample. All sources in the first 
< z < 0.1 bin have a redshift, and there is a small excess 
of pointed clusters between z — 0.8 and z = 1. In total, 188 
clusters out of the 347 ones selected for the CR-HR analy- 
sis have a redshift (66 of them being flagged as 'tentative' 
in our database, see Sect.[X]). However, including this par- 
tial redshift information in the MCMC analysis requires a 
precise knowledge of the associated selection process which 
is currently out of reach. 

6.2 Cluster X— ray profiles 

To account for the physical extent of clusters entering the 
survey selection function, we have introduced the x c $ pa- 
rameter linking the X-ray cluster extent to the cluster 
mass (x c ,o = r c /R 5 ooc)- Contrary to other studies, we 
did not assume a fixed physical size f or the core radius 
llPacaud et al. Il2007l; ISahlen et al.l l2009h nor a size distri- 
bution ( Burenin et alj I2007T ). In all configurations we in- 
vestigated, this parameter is constrained at the 10 — 20% 
level and is found to have a value of 0.24 ± 0.04 in the 
best-fit model. For a cluster of mass M200C = 10 14 /i _1 Mq 
at redshift z = 0.4, i?sooc is ~ 0.6 Mpc, and our result 
suggests a physical core radius of about 150 kpc, indeed 
typical of those found in other cluster studies at a simi- 
lar redshift (e.g. iPacaud et al. I l2007l ; IVikhlinin et al. 1998b; 
iBurenin et al.l l2007). As shown on Fig. 1191 there is a reason- 
able agreement between a /3-profile co mputed with our valu e 
of x Ci o and the AB-model used in (|Piffaretti et al.ll201ll) . 
base d on the local gas den sity profiles of the REXCESS clus- 



10.000 



1.000 



0.100 



0.010 



0.001 



AB model 
/S model 




base d on t he local g as den 
ters l|Croston et alj|2008h 





Figure 19. Comparison between th e AB model from 
iCroston et alj|2008t IPiffaretti et alj|201ll) and the /3-model used 
in this analysis (/3 = 2/3). These curves show surface brightness 
profiles normalized to the same total flux. The red curve corre- 
sponds to a /3 profile computed with the best-fit parameter x c ,o 
from the MCMC analysis. 



6.3 Scaling-laws evolution 

We parametrised the redshift evolution of each scaling law 
by two factors of the form (1 + z) 1 such that 7 7^ in- 
dicates a departure from self-similar evolution (Equ. [2] and 
[3j|. The modelling of the CR-HR distribution includes self- 
consistently selection effects and the evolution of cluster 
scaling relations which is a key point of such analyses (e . g . 
IPacaud et ai1l2007l ; IShort et al.|[2O10l ; iMantz et al.ll2010al ). 

We illustrate in Fig. the net effect of the combined 
factors E(z) and (1 + z) 1 in the evolution of scaling laws, 
using the best-fit values from our analysis. In this figure, the 
T-M relation has been computed by inverting Eq. [5] and is 
expressed in terms of the crticial mass M200C • The upper and 
lower boundaries have been computed using the covariance 
matrix of f2 m (which enters E(z)), 7 z ,mt and 7 z ,lt output 
of our MCMC chain. 

Our fit to the X-CLASS CR-HR diagram indicates a 
quasi non-evolving T-M200C relation (equivalently, a nega- 
tive evolution relative to the self-similar expectation), mean- 
ing that a cluster with a given mass M200C shows approxi- 
mately the same temperature at all redshifts. Our data also 
indicates a negative evolution of the L-T relation, below 
the self-simila r expectation, a result that has been found in 
simulations of lKav et alj (|2007l ) where feedback by AGN and 
stars is included in cluster simulations, and in IShort et al.l 
(2010) for their preheating model of clust e r evol ution. These 
trends are also observed bv lEttori et alj l|2004l ) for the M- 
T and L-T relatio ns b ut are i n conflict with s tudie s from 
iKotov fc Vikhlininl (e.g. 120051 ); iBranchesi et~al] (e.g. 120071 ). 
However, the comparison between these studies is hampered 
by the different selection processes entering different cluster 
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samples, which can have a tre mendous effect on the derived 
evolution (|Pacaud et alj|2007f ). 

Finally, we note that our results do not firmly exclude 
self-similar evolution both in the M-T and the in L-T rela- 
tions. The relatively large uncertainties on 7z ,mt and 7z,lt 
are due to the absence of redshifts and to the degeneracies 
with cosmological parameters inherent to the CR-HR anal- 
ysis (Fig. n~5|) . We expect the forthcoming XXL and XCS 
surveys to provide tighter constraints on the evolution of 
scaling laws thanks to the inclusion of cluster redshifts in 
both analyses. 



6.4 CR-HR outliers 

Fig[l6]shows the good agreement between our best-fit model 
and the CR-HR distribution of clusters in the sample. How- 
ever, we note the presence of sources outside the bulk of 
the diagram. The top-left par10 (CR < 0.06, HR > 1.4) 
contains regular, weak clusters for which measurement er- 
rors are large, particularly on HR (Fig. 1131) . Such hardness 
ratios could be explained by, e.g. a high metallicity of the 
intra-cluster medium favoring emission from bright spectral 
lines. The bottom-right regiorQ (CR > 0.1, HR < 0.7) con- 
tains 9 clusters whose morphology appeared strongly peaked 
at the centre, hence indicating a significant AGN or cool- 
core contribution; the net effect is to increase our integrated 
CR and modify the HR, with respect to normal clusters. 
One outlier (xclass 1937) is a compact group (HCG 057, 
iHickson. Kindl. fc Aumanl 19891 ) presenting a complex X-ray 
emission. A complete model would enclose cluster spectral 
peculiarities in the calculation of the CR-HR diagram; how- 
ever, we consider here that their number is sufficiently low 
to neglect their impact onto our results. 



6.5 Predictions for eRosita 

In paper I, we presented Fisher forecasts for an XMM 
100 sq. deg. cluster survey at a 10 ks depth, providing a sam- 
ple of 570 clusters. The CR-HR method was predicted to 
yield the following accuracy: cr(fi m ) = 0.09, a(ag) = 0.14, 
ct(7z,mt) = 0.6, <t( 7z , L t) = 2.3 and a(x C}0 ) = 0.04. The 
present study has been conducted for similar conditions 
(comparable area, a somewhat more stringent selection func- 
tion and partially deeper exposures). The good agreement 
between the predicted uncertainties and the output of the 
MCMC runs shows that both analyses are consistent with 
each other. We note however that there is not a one-to- 
one correspondence between them as the Fisher analysis in- 
volved more free parameters (though constrained with strin- 
gent priors) and there are 347 clusters in the present study, 
which are on average better measured than assumed in pa- 
per I. The difference in the total number of clusters is ex- 
plained by: i) the more stringent selection function, ii) the 
actual 90 sq.deg. coverage of the X-CLASS survey, iii) inho- 
mogeneities in the survey depth in terms of pointing back- 
ground and hydrogen column density and iv) differences in 



3 X-CLASS tags: 35, 86, 915, 997, 1032, 1655, 1741, 1886, 1947. 

4 X-CLASS tags: 102, 238, 541, 1020, 1218, 1480, 1906, 1937, 
2046, 2048, 2162, 2321. 



Table 8. Expected marginalized constraints on dark energy pa- 
rameters and parameters describing the scaling-law evolution for 
the 20,000 sq.deg. eRosita survey. In each case, the scaling rela- 
tions are let free in the analysis (slope, normalization and scatter). 
Planck priors are applied to the five parameters O m , <rg, fl^,, n B 
and h. 







CR-HR 


z- 


CR-HR 




No prior 


Planck priors 


No prior 


Planck priors 


wo 


0.6 


0.4 


0.1 


0.1 


W a 


1.0 


0.9 


0.3 


0.3 


7z,MT 


1.3 


0.1 


0.2 


0.05 


7z,LT 


0.8 


0.5 


0.3 


0.1 



the fiducial parameters of the Fisher matrix and those de- 
rived in the present work (e.g. x c fi — 0.1 in paper I and 0.2 
in this paper). 

The overall good agreement between predicted and 
measured uncertainties allows us to propose general p redic- 
tions for the Rosita all-sky survey l|Predehl et al.l l2010) to be 
obtained by the CR-HR method, following the formalism de- 
veloped in paper I. We assume a total area of 20,000 sq.deg. 
(extragalactic survey) and a custom selection function being 
a scaled-up version of the Cl + selection (Fig. [5]). Within our 
fiducial model, the survey is expected to yield 2.5 clusters 
per sq.deg., hence a total of 50,000 detected clusters. In con- 
trast to paper I, we do not assume priors on local scaling 
laws. We allow their normalization, slope and intrinsic scat- 
ter to vary, and let them evolve with redshift as (1 + z) 7 . We 
assume a prior Fisher matrix on fi m , ag, fit,, n s and h as 
will be available fro m Planck and calculated identically as in 
iPierre et al.l (|201ll ) (based on the Planck mission definition 
- BluebooiQ). 

We estimate measurement errors by assuming a mean 
exposure time of 2.5 ks and an effective area equal to that 
of XMM MOS+PN l|Predehl et al.ll2010l ). thus applying a 
factor 2 to the uncertainties quoted in paper I for a 10 ks 
XMM survey. We considered two extreme situations: either 
no redshifts are available and we apply a simple CR-HR 
analysis, or all clusters do benefit from photometric redshifts 
and we use the more fruitful z-CR-HR analysis with bins of 
Az = 0.03. Corresponding results are shown on Fig. I21l and 
quoted in Table [SJ for the dark energy parameters (wo,w a ) 
and for the parameters governing the scaling-law evolution 
(7z,mt, 7z,lt). We find that even without redshifts, the CR- 
HR method yields good constraints on the evolution of scal- 
ing laws, provided that the cosmological model is known at 
the accuracy expected from the Planck mission. The con- 
straints on wo and w a are less informative, but could be 
enhanced by a joint study of the angular correlation func- 
tion of the detected clusters. In particular, we notice that 
adding Planck priors to the analysis has a mild impact on 
the dark energy constraints, possibly because of the degen- 
eracies within the scaling-relation parameters and the rel- 
atively high uncertainties on the count-rate measurements. 
Adding redshift information considerably improves the ex- 
pected constraints on both sets of parameters. 



5 http:/ /www.rssd.esa.int/index.php?project=Planck 
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Figure 20. Scaling-law evolution as predicted by our model and parametrised by equations [2] and [3] The solid red line shows the best-fit 
model (Qm = 0.24, 7 Z mt = 0.8, 7 z ,lt = —1-3) with the dotted line indicating the 68% confidence boundaries. The black dashed line 
shows the self-similar expectation (7 Zj lt = 7z.MT = 0). 




Figure 21. Left: constraints on the dark energy parameters wo and w a for the eRosita 20000 sq.deg survey as predicted by the CR- 
HR method (blue), possibly supplemented by photometric redshifts for all clusters (z-CR-HR, green). No assumption has been made 
on scaling laws nor on their evolution and Planck priors were applied to O m , erg, f2 b , n s and h. The right panel shows the predicted 
constraints for cluster evolution 



7 CONCLUSIONS AND PERSPECTIVES 

We have presented the cosmological study of a sample of 
347 clusters detected in the full XMM-Newton archive us- 
ing X-ray criteria only (ancillary optical images were used 
to discard 'extended' sources irrelevant for our analysis such 
as nearby galaxies, saturated point-sources etc.). The study 
relies on the sole instrumental count-rates measured for each 
cluster in three X-ray bands. The selection function of our 
sample has been thoroughly defined by means of extensive 
image simulations and we proposed a method to account 
for the presence of pointed clusters in the analysis. We have 
then modeled the sample CR-HR diagram (whose proper- 
ties are discussed in paper I) by self-consistently including 



a ACDM cosmological model, X-ray scaling laws, selection 
effects and measurement errors. This allowed us to fit f2 m 
and as along with the parametrized evolution of scaling laws 
plus a parameter x c ,o characterizing the X-ray gas extent in 
clusters. We summarize below our main conclusions: 

• When setting the cosmological parameters to their 
WMAP-5 values, we observe a preference for the 'Non 
C ool Core' normaliz ation and slope of the L-T relation 
of iPratt et ail l|2009l ). if we assume an intrinsic scatter of 
a inh\T ~ 0.7. Our data then indicates a positive evolution 
of the M-T relation and a negative evolution of the L-T 
relation with respect to the self-similar expectation. 

• Fitting fi m and as together with the evolution of scaling 
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laws, we still find a preference for the 'NCC scaling law and 
find Q m = 0.24ijj;j$ and cr 8 = 0.88tg;i°, in agreement with 
the most recent cosmological studies. Again, the T-M and 
L-T scaling relations are found to evolve negatively with 
respect to the self-similar expectation. 

• Assuming the 'ALL' scaling relation, our data indicates 
fl m ~ 0.15 and as ~ 0.96 with a milder evolution of scaling 
laws. 

• The ad hoc parameter x C] o = r c /Rsooc giving the X- 
ray extent as a function of cluster mass is found to be well 
constrained within our framework, with a best value of 0.24, 
compatible with individual cluster studies. 

• The interpretation and use of the currently available 
local cluster scaling relations has proven one of the main 
hurdle of our study. Given that these relations do not agree 
with each other, it is probable that they have not been en- 
tirely corrected from the selection biases affecting the sample 
from which they are derived. 

• The scatter in the scaling relations plays an important 
role in the cosmological analysis and is probably degenerate 
with the slope and normalisation of the relations. It is likely 
that the scaling relations will be reliably determined only 
with very large cluster samples along with the simultane- 
ous fit of cosmology and selection effects. We stress that the 
CR-HR method is able to achieve this in a self-consistent 
manner, by-passing the tedious step that consists in deter- 
mining individual cluster masses. 

• As a logical follow-up of the present study, we propose 
predictions for the eRosita All-Sky survey. Assuming the 
Planck priors and letting all cluster scaling-law parameters 
free, we show that the z-CR-HR method will allow a de- 
termination of the equation of state of the dark energy a t 
the level of stage IV of the DETF l|Albrecht et al.| [2006). 
In addition, the cluster scaling-law evolution will be well 
determined. 

• The X-CLASS serendipitous cluster catalogue ex- 
tracted from the XMM archival data is available online at: 
http://xmm-lss.in2p3.fr:8080/14sdb71 
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APPENDIX A: X-CLASS CATALOGUE 
Al Database 

The X-CLASS catalogue is accessible through a dedicated 
database at http://xmm-lss.m2p3. fr:8080/14sdb/ 

It contains 845 CI clusters retained after data screening, 
in particular the 347 clusters used in the cosmological anal- 
ysis. The public part of the database contains 422 clusters 



selected identically as the cosmological sample (Cl + , with 
high- and low-cuts in [0.5-2] keV count-rate and hardness- 
ratio, see Sect. [3]) but extended up to 13 arcmin off-axis dis- 
tance (instead of 10 arcmin). Redshifts wer e obtained from 
the N ED and from the recent publication bv lMehrtens et all 
|201ll ) for the XCS survey. 

We describe below the informations contained in the 
database: 

• Object name: each cluster is referenced with a unique 
identifier ("tag"), a full name in the XMM-Newton format 
(XMMUJ) and a name output of the X-ray pipeline. The 
cluster name encloses the corresponding XMM Obsld (ex: 
0502430101) and truncated exposure time (lOks or 20ks) at 
which it has been detected. 

• Object position: the position of the cluster (right as- 
cension and declination) as provided by the pipeline is given 
in addition to the position measured by hand in the course 
of the cluster count-rate measurement (Sect. [2~4l . 

• NED identifications: column 'NED' lists all sources 
(galaxies, galaxy clusters, groups, QSO, etc.) within 3 ar- 
cmin of the cluster centre having a redshift indication from 
the NED (photometric or spectroscopic). An illustrative ex- 
ample is given in Table [ATI 

• Redshift: when a redshift indication is available, we 
provide a flag describing the current status of the redshift 
determination: 'confirmed' for a cluster definitely confirmed, 
'tentative' if less than 3 concordant redshifts within 3 arcmin 
are available and 'photometric'. Redshifts are quoted from 
NED first, then from the XCS-DR1 taking into account the 
provided flags. 

• X-ray properties: basic X-ray properties output of 
the XAmin pipeline (Sect. 12.21) are given for each cluster, 
in the [0.5-2] keV detection band: number of counts, total 
count-rate, apparent Extent and Extent Likelihood as well 
as the distance to the centre of the pointing it belongs to 
(off-axis). 

• Count-rates: count-rates measured manually in sev- 
eral energy bands are also available, in particular for the 
[0.5-2] keV band. In any case, count-rates are specified 'on- 
axis', i.e. corrected from the local exposure map, and do 
not include the filter and aperture corrections as discussed 
in Sect. 12.41 A webpage shows for each cluster the profiles 
generated for the count-rate measurement (see Fig. IA1I for 
a particular example). 

• Cluster images: X-ray photon images, filtered images 
and optical cut -outs from the Pa lomar Observatory Sky Sur- 
vey (POSS-II, iReid et"al1 Il99ll ) have been produced and 
linked to each catalogue entry. A screenshot image is shown 
on Fig. |H 

• Database interface: the electronic catalogue can be 
sorted according to any of entry of the database and can be 
downloaded as a machine-readable file. 



A2 Comparison to the XCS survey 

We compared our catalogue to the firs t release of the XCS 
survey XCS-DR1 |Mehrtens et al.ll201ll ) based on the analy- 
sis of all publicly available data in the XMM archive. For this 
comparison, we included our 347 clusters selected for the cos- 
mological analysis (with off-axis distance below 10 arcmin) 
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Cluster candidate : 0201090401 17 v3.3 cl lOks id : 908 
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Figure A2. Screenshot of a webpage from the X-CLASS database (http://xmm-lss.in2p3.fr:8080/14sdb/ ), for a particular cluster (tag 
908). Displayed arc the X-ray photon image (top-left) with associated contours, the wavelet-filtered (top-right) and gaussian-filtered 
(bottom-left) images and the optical (POSS-II) overlay (bottom-right) . 
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Table Al. Sample table attached to the 'NED' column in the database, for a particular cluster (tag 1758). This table lists all objects 
from the Nasa Extragalactic Database (April 2011) within 3 arcmin of the cluster center, with an associated redshift. Values in the 
table are directly copied from the NED and originate from various surveys/follow-ups. In particular, the redshift accuracy is highly 
inhomogeneous and in some cases a flag indicates its reliability. The velocity of each object is expressed in km/s. The distance is quoted 
relative to the cluster centre (units arcmin). 
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plus 75 clusters selected upon identical criteria, but with off- 
axis distances between 10 and 13 arcmin (thus 422 clusters 
in total). The XCS-DR1 sample is composed of 503 clusters 
which are optically confirmed and detected with more than 
300 counts in the [0.5-2] keV band. Clusters identified as tar- 
gets of a particular XMM observation are not included in the 
XCS-DR1. Fig. IA3I illustrates the comparison between the 
two samples. Out of our 422 clusters, 159 are new discover- 
ies (i.e. not in the XCS-DR1 and without information from 
NED). The differences between the two catalogues can be at- 
tributed to the different pointing selection and to differences 
in the X-ray detection algorithms. A further comparison be- 
tween the two analyses will provide useful insights into the 
different selection effects and possible systematics contained 
in both samples. 



APPENDIX B: BIAS CORRECTION FOR 
POINTED OBSERVATIONS 

We detail in this Appendix our method for correcting from 
the presence of pointed clusters in the XMM archive. As 
shown on Figures [10] and 1111 more clusters are detected 
in the inner 5 arcmin than expected from solely consider- 
ing the sensitivity gradient on the detectors. This sensitiv- 
ity difference is partly due to the EPIC vignetting func- 
tion (loss of ~ 60% in effective area at 10 arcmin off-axis 
("XMM-Newton Users Handbook", Issue 2.9, 2011 (ESA: 
XMM-Newton SOC)) compared to the centre). The degra- 
dation of the telescope PSF at larger off-axis distances is 
the second most important cause of sensitivity variation as 
it dilutes the signal of faint sources and distorts their shape. 



Bl Bias model 

We want to correct the CR-HR distribution from the bias 
due to pointed clusters, directly on the predicted CR-HR 



distribution. We assume that this correction does not de- 
pend on HR and divide the count-rate distribution in several 
bins indexed by j. We call rij the underlying cluster surface 
density, i.e. the value obtained after integrating the true 
cluster log(A r )-log(CR) in bin j. We introduce the survey 
selection function through the factor ej (comprised between 
and 1) such that the net number of clusters detected in 
bin j for an arbitrary region covering A = Q, f s k y sr. on sky 
is: 

Nf i =At i n i (Bl) 

Here the sky area is O = 47r(l — cos 70°) as we exclude the 
galactic plane (±20°) from the analysis. 

We now divide the survey in two sub-surveys labelled 
'in' and 'out'. The first one gathers all regions belonging 
to the inner [0-5] arcmin and the latter corresponds to the 
[5-10] arcmin regions of the same pointings. The total area 
covered by the survey reads A — A m + A out . Introducing /i™ 
(^ 1) the sensitivity of the inner sub-survey relative to the 
total survey we write (still for an unbiased sample): 

Ai-in * in in 

Nj = A ej n d rij 

We now consider that an unknown fraction fj of all 
existing clusters on sky has been observed and pointed in 
the central region of the detectors, thus augmenting AT™ by 
a quantity: 

Arin. pointed r~» j* in 

Nj 1 =Qrij fjtj (j,j 

Because those clusters have been removed from the sky 
population of clusters, the remaining density on sky is n'j = 
( 1 — fj ) rij and the survey provides a total number of clusters : 

a rtot .biased Arin, pointed . * I 

Nj = Nj v + Atj rij 

= ^ rij fj £j n? + A tj (1 — fj ) rij 

= ( 7 ^/ J Mi n + (l-./;)W t (B2) 

\Jsky / 
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Figure A3. Comparison between the X-CLASS and the XCS-DRl (Mchrtcns ct al. 2011) catalogues. Blue crosses stand for the cluster 
sample used in the cosmological analysis, blue squares for the "extended" cosmological sample (i.e. up to 13 arcmin off-axis distance on 
the XMM detectors). Green diamonds show the 92 sources of our sample located at the centre of the XMM field of view (i.e. less than 3 
arcmin off-axis distance). 



Comparing equations IB 1 1 and IB2I we obtain the bias 
factor: 



tot, biased /7\7-tot ^ 

fsky 



/Nf* = -^fitf + (l-fi) (B3) 



B2 Bias estimation 

Following previous equations, the expected number of clus- 
ters in the 'in' survey writes: 



r-\ £ in i a in in 

= " rij fj e 3 //, + A t j ii j 



1 A in 

_Jsky -rY 



(B4) 



while the expected number of clusters in the 'out' survey is: 

Nj 0t (B5) 



jyout, biased ^ y ^ 



, in A 



At this point, only two quantities are unknown, fj and 
Nj 0t . Parameters A and A ln directly come from the survey 
geometric design. The factor ptf is obtained by comparing 
the results of the simulations (Sect. [3]) for clusters in the 



full [0-10] arcmin off-axis area and for clusters in the central 
[0-5] arcmin region. 

We compute independently in each bin the joint likeli- 

out - — ~in 

hood for the observed (biased) quantities Nj and Nj on 
a fine bidimensional grid sampling values for fj and Nj° . 
We then compute the marginalized probability distribution 
P(fj) assuming flat priors for JVj ot by numerical integration 
of the sampled likelihood. 

The expectation value and variance for Fj are finally de- 
rived by integrating Eq. IB3 1 against P(fj) and are displayed 
on Figure Q2] for the present sample. The predicted, unbi- 
ased CR-HR distribution is multiplied by the expectation 
value of Fj linearly interpolated at each CR value. 

This model relies on the hypothesis that all pointed 
clusters are centered onto the detectors. It thus neglects 
spatial correlation effects which may artificially boost the 
number of clusters in the surroundings of pointed clusters, 
in particular in the outer parts of the detectors. 



This paper has been typeset from a TeX/ ETfrjX file prepared 
by the author. 
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Figure Al. Screenshot of a webpage from the X-CLASS 
database (http://xmm-lss.in2p3.fr:8080/14sdb/ ), for a particular 
cluster (tag 2094). Displayed are the cumulative, background- 
subtracted, count-rate profile (top) the differential count-rate 
profile (2nd panel), the corresponding signal-to-noise curve (3rd 
panel) and the surface brightness profile (bottom). The blue verti- 
cal lines shows the radius in which the measurement is performed. 
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